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ABSTRACT 

We obtained a time series of spectra covering the secondary maximum in the /-band of 
the bright Type la supernova 2014J in M82 with the TIGRE telescope. Comparing the 
observations with theoretical models calculated with the time dependent extension of 
the PHOENIX code, we identify the feature that causes the secondary maximum in the 
/-band light curve. Fe II 3/ 6 ( 3 /3)4s —3/ 6 ( 5 /))4p and similar high excitation transitions 
produce a blended feature at ~ 7500 A, which causes the rise of the light curve towards 
the secondary maximum. The series of observed spectra of SN 2014J and archival data 
of SN 201 lfe confirm this conclusion. We further studied the plateau phase of the R- 
band light curve of SN 2014J and searched for features which contribute to the flux. 
The theoretical models do not clearly indicate a new feature that may cause the 17- 
band plateau phase. However, Co II features in the range of 6500 — 7000 A and the 
Fe II feature of the /-band are clearly seen in the theoretical spectra, but do not 
appear to provide all of the flux necessary for the /7-band plateau. 

Key words: radiative transfer - supernovae: general - supernova: individual: 
SN 2014J. 


1 INTRODUCTION 

The light curves of Type la supernovae (SNe la) are of gen¬ 
eral interest, in particular for applications in cosmology. The 
optical li ght curves of SNe la admit a lightcurve shape re¬ 
lation ([Phillipsl 19931 : iRiess et all 1 19961 : IPhillips et alJll999t 
iGoldhaber et al] 200lh which allows them to be used as 
“correctable candles” and hence, their distances can be de¬ 
termined from the observed photometry. The underlying 
physical explanation for the light curve shape relation is 
that the total 56 Ni mass produced in the explosion de¬ 
termines to first order the temperature and thus, more 
56 Ni leads to higher temperatures, wh ich leads to higher 
opacities, increasing the diffusion time ([ Nugent et al. 1995|: 
KhokhlovetaLl 19931: Hoflichet U1 jR95lJPintofcEastmanl 


2000al lbl: Pinto fe Eastman 2001; Mazzali et al. 2 00Hl . It has 


been questioned whether the primary effect of the 56 Ni mass 
is on the diffus ion time or rather on the color evolution 
of the spectra (K asen fe Wooslevl l2007l l. but the ultimate 
source for the observed luminosity of a SN la is the decay 
of 56 Ni, and thus to zero order the luminosity should be 
correlated with the mass of 56 Ni. Using the light curves of 
SNe la for distance measurements, led to the discovery of 
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the phenomenon that the Universe’s exp ansion is acceler¬ 
ated, i.e, the discovery of the dark energy llRiess et al.lll998l : 
IPerlmutter et all 1 1 9991 ) . The nature of the dark energy is of 
fundamental importance for cosmology. 

Many physical details of Type la supernovae events are 
still not fully under stood I'see lHillebrandt fc NiemeveifcOOOl ; 
IParrent et al.ll2014i . for reviews). One of the most important 
open questions is the nature of the progenitor(s) of SN la ex¬ 
plosions. The progenitor is believed to be a CO white dwarf 
in a binary system. Its companion star could be a main se¬ 
quence star, an evolved star, or another white dwarf. The 
primary white dwarf grows in mass via accretion, which can 
be slow in the Chandrasekhar mass scenario, or rapid in 
scenarios that involve the merging of two white dwarfs, (see 
iLevanon et al .1120151 : iMaoz et al .11 20141 . for reviews of the sug¬ 
gested possibilities). In fact, SNe la may come from a variety 
of progenitor systems, and detailed observations of nearby 
SNe la may lead to the understanding of the nature of thes e 
progenitor systems llNueent et al.ll201ll : lDildav et al.l 1 2 01211 . 
There is further ongoing discussion about the explosion 
mechanism of Type la supernovae. An instantaneous deto¬ 
nation in Chandrasekhar-mass WDs has been ruled out, be¬ 
cause it fails to produce intermediate mass elem ents (IME), 
which are observed in the spectra of SNe la ( Arne ttl ll969h . 
In sub-Chandrasekhar-mass WDs an instantaneous detona- 
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tion can produce IME due to the low er central densities 
l|Shigevama et alii 19921 : ISim et alX lOlO). However, the mass 
of the WD is rather tightly constrained to produce both 
iron group elements and IME and the low densities will 
not produce observed stable iron group elements. Several 
different explosion models such as a pure deflagration, de¬ 
layed detonations, gravitationally confined detonations, vio¬ 


lent mergers J _and collis i ons have bee n sugg ested jNomotc 
1984; iKhokhlovLJlMjyPlewartalJ J^OOJ^Jjordanet^l. 


II 

2008 L 2012^. Pakmoret_ah[ 2012. 2013; iRosswog et al. 20091 : 
Kushnir et al.l 20131 : Long et al.l 12014 1. which could explain 


the appearance of the intermediate mass elements observed 
in SNe la spectra. To understand the physical processes in 
a Type la supernova further, different numerical codes have 
been developed and used to model the spectra and light 


Maeda et al. 2014 

: Wollaeeer & van Rossum 2014; Hoflich 

2003 

; Kasen et al. 

2006: Sim et al. 20llS: Dessart & Hillier 


We focus here on the light curves in the near-infrared 
wavelength region. The light curves of the near-infrared 
bands 7, J, 77, and K very often show a secondary max¬ 
imum (lElias et al.lll98ll : lHoflich et al.ll2002f) . In his detailed 
study of near-infrared light curves of SNe la, iKase 3 (|2006h 
showed that the secondary maxima are caused by the recom¬ 
bination of iron peak elements, which comes from a drastic 
opacity change between ionization stages III to II of the iron 
elements. This observation has been confirme d wi th other 
theoretical radiative transfer calcul a tions lljack et al. 


iDessart et al.ll20l3 : iGall et alll2012l ) . Ijack et al.l 1 2012 ) were 


120121 : 


able to assign specific features of Fe II and Co II that cause 
the individual secondary maxima in the 7, J, 77 and K 
bands. However, their investigatio n was based only on the¬ 
oretical light curve calculations. iDessart et al.l (l2014l f found 
that [Co II] lines play a strong role, but this result has not 
been confirmed in other calculations and may be due to the 
way that levels and superlevels are coupled in their approx¬ 
imations. 

In this work, we will present observational evidence that 
an Fe II feature, produced by a blend of Fe II lines, causes 
the secondary maximum in the 7-band and additionally, we 
present an investigation of the plateau phase of the 7?-band 
light curve. In §[2j we present the details of the observations 
of SN 2014J and describe the methods used to calculate the 
theoretical light curves with the PHOENIX code. In § [HJ we 
present the results obtained for the secondary maximum in 
the 7-band and the plateau phase of the 7?-band. 


2 METHODS 

We obtained a time series of observations of high resolu¬ 
tion spectra of the bright SN la, 2014J in M82 with the 
TIGRE telescope. We fitted the model light curves calcu¬ 
lated with the PHOENIX code to the observed light curves in 
the near-infrared 7 and 7?-bands. The observed and theo¬ 
retical spectra are then compared to determine the features 
that cause the secondary maximum in the 7-band and the 
plateau phase in the R-band. 

SN 2014J was discovered bv lFossev et al.l (l2014l 'l in the 
nearby galaxy M82 approximately 1 week after explosion. 


We use f or the explos ion time January 14.75 UT as found 
bv IZheng et al.l (120141 '). 


2.1 Observations 

We observed SN 2014J in M82 with the TIGRE telescope 
situated in Central Mexico close to the city of Guanajuato 
dSchmitt et al .lEfOlll ). The telescope TIGRE is a 1 .2 m mir¬ 
ror telescope that operates completely robotically. Its orig¬ 
inal design was to monitor the stellar activity of solar-like 
stars, but other observational campaigns are also feasible. 
The TIGRE telescope is equipped with the HER.OS spectro¬ 
graph, which has a high resolution R ~ 20,000 and covers 
the wavelength range from about 3800 A to about 8800 A. 
The spectrum is divided into two channels (red and blue). 
The data reduction pipeline is also fully automatic. 

Beginning with the first observation during the night of 
23rd of January 2014, we were able to observe a spectrum 
of SN 2014J i n M82 almost every night following the discov¬ 
ery. See [jack et al.l (l2015l l for a complete presentation of the 
time series of SN 2014J spectra observed with the TIGRE 
telescope. We used an exposure time of 3 hours for every 
observed spectrum. The spectra were obtained with the full 
spectral resolution R ~ 20, 000. Since a supernova spectrum 
shows only very broad features and such a high resolution 
is not necessarily required, we binned the spectrum down to 
a resolution of 10 A. The binning significantly improves the 
signal to noise ratio. For this work, we are interested in the 
near-infrared wavelength range, so that we will only use the 
spectra of the red channel (5800 A to 8800 A). SN 2014J 
showed significant reddening so that the obtained spectra in 
the blue channel (3800 A to 5800 A) observed at such high 
resolution have a much lower signal to noise. 


2.2 Theoretical light curves 

In order to model the theoretical light curves in different 
bands of SN 2014J, we used the time dependen t extension 
of the PHOE NIX code as described in detail in Ijack et ahl 
J2009I . l201lf) . The code follows the evolution of the super¬ 
novae envelope after the explosion enforcing energy conser¬ 
vation. This includes energy deposition due to y-rays from 
the radioactive decay of 56 Ni and 56 Co. Cooling due to the 
adiabatic expansion of the envelope, which is assumed to be 
homologous, is also included. For the transport of energy 
by radiation throughout the envelope, we solve the spheri¬ 
cally symmetric special rela tivistic radiative tran sfer equa¬ 
tion (lHauschilddll99l3 : lHauschildt &: Baronl fl999h . The en¬ 
velope was divided into 128 layers. About 1,000 time steps 
are performed for each point in the light curve. We usually 
calculate a point in the light curve every 1 or 2 days. All 
light curves have been calculated assuming local thermody¬ 
namic equilibrium (LTE) and using about 2000 wavelength 
points. The final spectra presented in this work have been 
calculated using a higher resolution of about 20,000 wave¬ 
length points. PHOENIX uses the update Kurucz database Q 
which contains about 80 million atomic line transition. For 
the initial input structure, we used the results of the W7 


1 http://kurucz.harvard.edu/linelists/gfall/ 
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deflagration model (lNomotolll984h . which also includes the 
specific non-homogeneous abundances. 

In this work, we focus on the secondary maxima in the 
near-infrared wavelength region. We seek to obtain the best 
fit to the individual observed light curves of SN2014J in 
the respective bands. As described in ljack et al.l (|2012| ), with 
the above assumptions, it is necessary to vary the equivalent 
two-level atom thermalization parameter when solving the 
radiative transfer problem. The source function of the radia¬ 
tive transfer equation including scattering for an equivalent 
two level atom can be written as 

= (1 — ex) J\ + e\B\. (1) 

5a is the source function, B\ is the Planck function, and 
Ja is the mean intensity. All these quantities are wavelength 
dependent. For the LTE PHOENIX calculations, it is possible 
to set a wavelength independent factor t = t\ — constant to 
approximate LTE line scattering over the whole wavelength 
range. We use this method here for computational expedi¬ 
ency. Since the SN la envelope becomes thinner during its 
free expansion phase, scattering becomes more and more 
important, and the thermalization parameter e decreases. 
Applying this method of a decreasing thermalization pa¬ 
rameter, we obtained fits to the observed light curves of SN 
2014J. 


3 SECONDARY MAXIMUM 

We compare the spectra of SN 2014J observed with the TI- 
GRE telescope with the theoretical spectra calculated with 
PHOENIX focusing on the secondary maximum of the /-band 
light curve and the plateau phase observed in the /?-band. 
All observed spectra and light curves have been dereddened 
with the values E(B — V) = 1.33 and Ry = 1.3 found by 
Amanullah et al. (2014) in their study of the extinction law 
of SN 2014J. 

3.1 /-band 

We obtained a time series of spectra of SN 2014J during the 
rise of the light curve towards the secondary maximum in 
the /-band. Unfortunately, the HEROS spectrograph does 
not cover the whole range of the /-band filter, meaning that 
we could not reconstruct a light curve from our observed 
spectra. Therefore, we used the observed light curve of the 
AAVSO database^ to find the best fit of the model light 
curve calculated with PHOENIX. As described above, we var¬ 
ied the thermalization parameter e to obtain the best fit to 
the observed light curve in the /-band. Although this is not 
a physically correct method to calculate a light curve, the 
overall trend of a decreasing e is physically motivated. How¬ 
ever, the goal here is to use this fit to study the spectral 
evolution of a SN la towards the secondary maximum. 

In Figure |T[ the filled circles represent the light curve 
of SN 2014J in the /-band obtained from AAVSO observa¬ 
tion data. The observed light curve has been dereddened. 
We also plotted the AAVSO light curve in absolute magni¬ 
tude using the distance of 3.5 Mpc to M82, which results in 

2 www.aavso.org 



Figure 1. The circles represent the dereddened AAVSO observed 
light curve in absolute magnitude of SN 2014J in the /-band. The 
solid line shows the best fit obtained from PHOENIX models. The 
values for e are indicated along the bottom of the plot and the 
vertical dotted lines demarcate the range. 



wavelength [A] 

Figure 2. Theoretical spectra at different epochs after the ex¬ 
plosion calculated with PHOENIX. At ~ 7500 A, a new feature of 
blended Fe II lines appears and causes the secondary maximum. 
The dotted line indicates the /-band filter response. 


a distance modulus of 27.7 mag, which is comparable with 
the value of 27.6 mag found bv lFolev et alJ ll2014l) . The solid 
line represents the best fit of the theoretical light curve cal¬ 
culated with PHOENIX models. As stated above, the models 
were calculated with different values of e for different epochs. 
This causes the kinks in the model light curve. The rise to¬ 
wards the secondary maximum is clearly visible and starts 
around 32 days after the explosion. The secondary maxi¬ 
mum is reached at around 45 days. These are the epochs 
that are interesting for studying the spectral evolution and 
identifying the feature that causes this secondary maximum. 

In Figure [2] we show a plot of the theoretical spectra 
calculated with PHOENIX at different epochs obtained from 
the best fit to the observed light curve. The spectra are 
shown in flux received at Earth if the models would be at 
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wavelength [A] wavelength [A] 


Figure 3. Spectra of SN 2014J observed on two different days 
before and during the secondary maximum. The appearance of 
the feature of blended Fe II lines can be seen in the observations 
at ~ 7500 A as well. The dotted line represents the I -band filter 
function. 


Figure 4. Spectra of SN 201 lfe observed. The appearance of the 
feature of blended Fe II lines can be seen in the observations at 
~ 7500 A as well. The dotted line represents the /-band filter 
function. 


the distance of SN 2011fe (6.6 Mpc) to be able to com¬ 
pare the fluxes directly. The dotted line shows the /-band 
filter response. At that epoch the spectrum does not show 
many clear features. The brightness decreases after the / 
maximum. The Ca II IR triplet feature at 8500 — 9000 A 
becomes prominent. However, this feature, does not show a 
change in brightness during the rise to the secondary maxi¬ 
mum. From the beginning of the rise (32 days) up to the 
secondary maximum (45 days) there appears a new fea¬ 
ture of Fe II at a wavelength of ~ 7500 A. In Fe II there 
are numerous transitions from relatively high excited states 
in this wavelength region, for example transitions from 
3d 6 ( 3 D)4s —3d 6 ( 5 D)4p, 3d 5 ( 4 P)4s4p( 3 P°)-3d 6 ( 5 //)6s, and 
similar transitions, which have Einstein coefficient A 21 val¬ 
ues > 1 x 10 7 s~ 4 . Thus, this feature causes the observed 
secondary maximum in the /-band of SNe la. 

The observed spectra also clearly show the appearance 
of the Fe II feature mentioned above. Figure [3] shows two 
observed spectra of SN 2014J in the same wavelength range 
as was used for the theoretical spectra in Figure [2] The dot¬ 
ted line represents the /-band filter response. The spectrum 
from the 16th of February was taken directly before the be¬ 
ginning of the rise towards the secondary maximum. The 
spectrum from the 1st of March was observed during the 
secondary maximum. These observed spectra correspond to 
the theoretical spectra of day 32 and day 45 in Fig. 2. The 
feature of Fe II at around 7500 A also appears in the ob¬ 
served spectra. Note that there is a strong telluric absorption 
feature around 7600 A, which contaminates the spectrum in 
this region. However, we can clearly confirm that this Fe II 
feature actually causes the secondary maximum in the I- 
band of SN 2014J, and thus, it is likely the cause in most 
SNe la. 

Since our observed spectra of SN 2014J do not cover the 
whole wavelength range of the /-band, we used archival data 
to obtain additional evidence that the Fe II feature causes 
the secondary maximum. 

In Figure [4] we show two spectra of SN 2011fe. The 


spectra were observ ed by the Nearby Supernova Factory 
ll Pereira et al ] ml) in a time series of spectra. The spec¬ 
trum from the 24th of September 2011 coincides with the 
minimum before the rise towards the secondary maximum in 
the /-band. The secondary maximum was observed around 
the 8th of October 2011. Unfortunately, there does not exist 
a spectrum of that day. Therefore, we used the spectrum 
of the 4th of October. In the plot one can clearly see that 
a feature around 7500 A arises. This gives more observa¬ 
tional evidence that this Fe II feature causes the secondary 
maximum in the /-band of SNe la. 

3.2 //-band 

The light curve of a typical type la supernovae in the R- 
band does not show a secondary maximum, but usually a 
short plateau phase. Since all other light curves in the opti¬ 
cal wavelength range show a steady decline after maximum 
light, there should also be a feature that causes the effec¬ 
tive //-band photosphere to “hang up”. We used the same 
method to identify this feature as in the study of the /-band 
secondary maximum. Using a varying line thermalization 
parameter e, we obtained a best fit to the observed (dered- 
dened) //-band light curve of SN 2014J. 

In Figure[S] the filled circles show the observed AAVSO 
light curve in the //-band. Between about 30 and 40 days 
after explosion the //-band light curve of SN 2014J shows a 
short plateau phase, which we also modeled with the the¬ 
oretical light curve of PHOENIX. The solid line represents 
the best fit obtained from PHOENIX models by varying the 
thermalization parameter. Note that the variations of e are 
slightly different from those used to compute the /-band 
light curve. The differences, while small, represent a limita¬ 
tion in our modeling approach. However, the overall trend 
of a decreasing e is the same for all near-infrared bands. 
Scattering becomes more and more important as the SN la 
envelope expands. 

In Figure [6] we compare the //-band light curves that 
have been calculated using the values for e used for the R- 
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Figure 5. The circles represent the dereddened AAVSO observed 
light curve in absolute magnitude of SN 2014J in the //-band. The 
solid line shows the best fit obtained from PHOENIX models. The 
values for e are indicated along the bottom of the plot and the 
vertical dotted lines demarcate the range. 


Figure 7. Spectra of the best fit to the //-band light curve. Two 
features of blended Co II lines and blended Fe II lines contribute 
to the flux that causes the //-band plateau phase. The dotted line 
represents the //-band filter function. The model luminosity has 
been scaled by the distance to SN 201 lfe. 



Figure 6. Here we show the PHOENIX //-band light curve calcu¬ 
lated with the e values for the /-band and //-band. 


band and for the /-band. One can see some clear differences 
during the phase of the minumum before the plateau phase. 
However, since e is wavelength dependent it is not possible to 
distinguish between a possible inaccuracy of the models or a 
real difference in the value of e for the different wavelength 
regions. 

We then studied the spectral evolution during this ob¬ 
served plateau phase. In Figure |T] we show the theoret¬ 
ical spectra calculated with PHOENIX at the beginning of 
the plateau phase (30 days) and at the end of the //-band 
plateau (40 days). The dotted line represents the //-band fil¬ 
ter function. After the R maximum (16 days), the brightness 
decreases in the whole wavelength range. At 30 days, which 
marks the beginning of the plateau phase, the characteristic 
P-Cygni Si II feature around 6000 A to 6500 A strength¬ 
ens. This feature changes between 30 and 40 days during 
the plateau phase. However, the contribution of this feature 
to the brightness in the //-band decreases and can, there¬ 



wavelength [A] 


Figure 8. Spectra of SN 2014J observed on two different days at 
the beginning and the end of the //-band plateau phase. The Fe II 
features is clearly visible while the Co II feature of the models has 
not been identified. The dotted line represents the //-band filter 
function. 


fore, not cause the plateau phase. In the models, blends of 
high excitation Co II lines appear in the wavelength region 
6500 — 7000 A, however the total flux in the models is not 
large enough to produce the observed plateau. The feature 
of Fe II at 7500 A that causes the secondary maximum in 
the /-band is also present in the //-band filter, but again 
the filter sensitivity is falling off here so the total flux is not 

large enough to produce the observed pl ateau. __ 

While nickel holes are seen in SNe la iFesen et alj2007h . 
the nickel hole in W7 is too large to be compatible with 
galactic nucleosynthesis and thus, the center of W7 has too 
much stable Ni and Fe, which will reduce the total amount 
of Co. However, the Co II lines seen in the models are not 
easily identified in the observations (see below). 
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Figure [8] shows the observed spectra of SN 2014J in the 
R-band using the same wavelength range as for the theoret¬ 
ical spectra of Fig. [T] The dotted line once again represents 
the R-band filter function. We plotted the spectrum from 
14th of February, which corresponds to the beginning of the 
plateau phase (30 days). The end of the plateau phase at 
40 days is represented by the observed spectrum from the 
24th of February. One can clearly see the change of the ex¬ 
pansion velocity of the characteristic Si II feature as the 
photosphere moves inwards. There is no clear feature aris¬ 
ing in the R-band that obviously causes the short plateau 
phase. Of course, the Fe II feature at around 7500 A also 
contributes somewhat to the light curve in the R-band. As 
discussed above Co II lines in the 6500 — 7000 A range are 
seen in the models, but their signatures are not strong in 
the observed spectra. 


4 CONCLUSIONS 

We observed a time-series of spectra of the nearby Type la 
supernova 2014J in M82 with the TIGRE telescope. This 
time-series also covered the secondary maximum in the I- 
band and the plateau phase in the R-band. We calculated 
the respective theoretical light curves with the time depen¬ 
dent extension of the PHOENIX code by varying the line ther- 
malization parameter to obtain a best fit to the observed 
light curves. For the /-band light curve of SN 2014J, we 
find that the secondary maximum is caused by a blend of 
Fe II high excitation lines because of the recombination of 
iron peak elements from ionization stage III to stage II. This 
confirms previous work studying near-infrared light curve s 
of SNe la (lKasenl2006l : IJack et al.f2012l : lDessart et al Jl 201 4lf . 
Studying the spectral evolution of SN 2014J, we find obser¬ 
vational evidence that the Fe II feature at around 7,500 A 
causes the secondary maximum in the /-band. While W7 
is only a parameterized model of a real SN la, the interior 
composition is fairly generic due to the nuclear physics and 
thus, our conclusion should be fairly general. However, the 
details are a bit more complicated, and recent explosion cal- 
culations show quite some diversity in the ejected materia l 
llFink et al.ll2014l : IWooslev fc Kasenll201ll : iMoll et afll/OUl) . 

As far as we know, there exists no detailed study of the 
R-band light curves of Type la supernovae. We observed 
SN 2014J, and it — like most other SNe la — shows a short 
plateau phase in the //-band between 30 and 40 days after 
the explosion. We reproduced the observed light curve with 
our PHOENIX models to be able to identify the feature that 
causes this short plateau phase. Co II lines in the range 
6500 - 7000 A do increase in flux in the models, and the 
Fe II lines to the red may also contribute somewhat, but 
the total flux from these features does not appear to be 
large enough in the models and the Co II signatures are not 
clearly seen in the observations. 

For future work, it will be very interesting to study the 
spectral evolution of a Type la supernova in the other near- 
infrared bands J, H, Y, and K. With this one might be able 
to identify the features which cause the second ary maxima 
in these bands as proposed by I Jack et al.l d2012tl . A further, 
more detailed study of the R-band is important, since under¬ 
standing the detailed composition and ionization structure 


of the ejecta with time will help to constrain explosion mod¬ 
els. 
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